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Emission  of  Particles  From  a 
Charged  Sphere  Into  a  Magnetic  Field 


1.  INTRODUCTION 

1-3 

A  number  of  experiments  dealing  with  the  emission  of  beams  of  charged 
particles  from  material  bodies  Into  the  earths  magnetosphere  have  been  reported 
in  the  literature.  Since  many  of  the  calculations  relating  to  these  experiments 
attempt  to  anticipate  and  to  more  or  less  model  actual  ambient  conditions,  they 
are  so  involved  that  little  insight  into  the  basic  mechanisms  that  influence  the 
emission  process  is  gained.  To  fill  this  gap,  a  prototype  problem  —  namely,  to 
calculate  in  the  presence  of  a  spherically  symmetric  electric  field,  the  effect  of 
an  assumed  constant  magnetic  field  on  particles  emitted  normally  from  a  sphere  — 
has  been  solved.  In  particular,  if  3  is  the  ratio  of  cyclotron  to  sphere  radius  and 
0  the  angle  between  the  direction  of  emission  and  that  of  the  magnetic  field,  rela¬ 
tionships  between  3  and  0  which  separate  all  emitted  particles  into  two  classes  - 
those  that  return  to  the  sphere  and  those  that  do  not  —  a;'e  obtained.  The  result  of 
this  exercise  is  a  set  of  curves,  one  for  each  sphere  potential,  that  give  insight 
into  the  role  of  the  magnetic  field  in  general  emissive  processes. 
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Calculations  of  actual  charge  densities  and  potential  distributions  surround¬ 
ing  the  sphere  require  self-consistent  solutions  of  charged  particle  trajectory  and 
Poisson  equations,  a  lengthy  and  difficult  procedure.  Here  we  do  not  carry  out 
these  calculations,  but  rather  obtain  qualitative  indications  of  the  effect  of  space 
charge  by  assuming  potentials  proportional  to  1  /r*  (  r'  =  radial  coordinate)  and 

determining  separation  curves  for  each  of  several  integral  values  of  N.  Although 

not  all  potentials,  $(r'),  can  be  represented  by  a  uniformly  convergent  power 
oo  a 

series,  L  — =r*  ,  in  the  interval  r  <  r'  <  «,  one  can,  by  summing  terms,  obtain 

n-0  r'n  ° 

information  on  a  fairly  large  class  of  potentials. 


Z  CALCULATIONS 


If  we  assume  a  potential  of  the  form 


♦  ■♦.of  • 


(1) 


the  appropriate  equations  and  initial  conditions  for  this  problem  are  the  following: 


__  d^  r*  /  .  r 1  r^  $  \ 

dt2  =  6  (hT  X®  *N  r|N+2°  ) 


r *  =  T  r  sin  0  +  I?  r  cos  0 
J  o  o 


(2) 


dr' 


-St-  =  Tv  sin  0  +  Tc  v  cos  9 
at  v  o  o 


t  =  0 


(3) 


Here 


m,  e  =  particle  mass  and  charge  respectively 
P  =  radius  vector 
r»  =  |p| 

fi  =  magnetic  field  vector 

rQ  =  sphere  radius 


=  sphere  potential 

T.Tt  *  unit  vectors  along  rectangular  y  and  z  axes  respectively 
0  =  angle  between  magnetic  field  and  emission  direction 

vQ  =  initial  velocity 
t  =  time  . 

By  the  following  changes  in  variables. 


r  =  r*/r 


eB 

_  O  4. 

T  =  — - t 

m 


m  v 

^  =  r-Fir 

o  o 


m*o 

a  2  Z2 
r  e  B 
o  o 


(Bq  =  |b|),  Eq,  (2)  may  be  put  into  the  rectangular  component,  non-dimensional 
form: 


.2  . 

d  x  dy  X7  x 

7J  3r  ‘  tfN  -NTS’ 

dr  r 


d_y  _  dx  M  y 
—h  -  -  SF  •  ffN  ~n+3 

or  r 


(4) 


d^Z  mj  2 

-«N  ^5 


with 


x  =  0 

y  =  sin  0 

z  =  cos  0 

T  =  0 

■g?  =  0  sin  6 

dz  « 

-p-  =  0  cos  0 

as  corresponding  initial  conditions. 
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If  these  equations  are  converted  to  an  equivalent  set  of  six  first-order  equa¬ 
tions,  they  may  be  solved  by  use  of  standard  computer  codes.  To  generate  the 
curves  that  separate  escaping  from  returning  trajectories,  values  of  a  and  0 
(subject  to  the  energy  constraint  0  =  J2a)  are  chosen.  Starting  with  9  &  v/2, 
which  guarantees  a  return,  trajectories  are  taken  for  values  of  9  successively 
decreased  by  0.05  radian  until  an  escape  is  obtained.  Starting  with  the  value  of 
9  at  which  an  escape  is  obtained,  trajectories  are  taken  for  values  of  B  successively 
increased  by  0.005  radian  until  a  return  is  obtained.  This  process  is  continued 
until  the  desired  accuracy  of  the  value  of  9  is  obtained. 


3.  RESULTS 


The  results  for  N  =  1  are  shown  in  Figure  1.  For  any  pair  of  values  (0,  0) 

falling  above  a  given  curve,  the  particle  escapes;  for  any  pair  falling  below  a 

given  curve,  the  particle  returns  to  the  sphere.  With  no  magnetic  field,  the 

separation  curves  would  be  continuations  of  the  straight  horizontal  lines  out  to 

9  -  n/2;  these  lines  are  representations  of  the  dimensionless  energy  relationship 
o 

a  =  0  /2.  The  portions  of  the  curves  that  deviate  from  the  straight  lines  thus  give 
the  effects  of  the  magnetic  field  in  returning  particles  to  the  sphere. 

For  a  =  0,  the  trajectory  is  readily  solved  for  analytically.  We  have,  for 
this  case, 

r2/r2  =  1  +  2  02  sin2  0  +  20  sin2  0  sin  t  -  2  02  sin2  0  cos  r 
+  02  r2  cos2  0  +  2  0  r  cos2  0 

2  2 

The  condition  that  particles  be  recollected  is  that  the  equation  r  /r  =  1  have  a 
solution  for  some  r  >  0.  The  relationship  between  0  and  0  that  separates  escaping 
from  recollected  particles  will  then  be  given  by  eliminating  t  from  the  simultaneous 
implicit  equations; 

F  =  2  tan2  0  (0(1  -  cos  t)  +  sin  t]  +0t2  +  2t  =0  (5) 

-j  -gp  =  tan2  0  (0  sin  r  +  cos  t]  +  0  t  +  1  =  0  ,  (6) 


where 


F  =  <r2/r>o  - 


1 

2  • 
cos  0 
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Figure  1.  3,  Ratio  of  Cyclotron  to  Sphere  Radius,  vs  0, 

Angle  Between  the  Direction  of  Emission  and  That  of  the 
Magnetic  Field,  for  Various  Values  of  a  and  N  =  1.  The 
curves  separate  all  emitted  particles  into  those  which  re¬ 
turn  to  the  sphere  (below  curve)  and  those  which  do  not 
(above  curve) 

For  large  values  of  3,  Eqs.  (5)  and  (6)  can  be  solved  analytically  to  obtain 
tan  0  =  2  ?r3  t 

which,  for  large  3,  agrees  well  with  numerical  solutions  of  Eqs.  (5)  and  (6)  and 
also  with  solutions  of  Eqs.  (4)  with  a  =  0.  It  is  seen  from  Figure  1  that  for 
a  =  0  there  is  an  angle  6  s  1. 13  radians  below  which  all  particles  escape  indepen¬ 
dent  of  initial  energy. 

The  simulated  effects  of  space  charge  are  demonstrated  in  Figures  2,  3,  and 
4  where,  for  values  of  a  =  0.01,  1.0,  and  100,  separation  curves  for  various 
values  of  N  are  plotted.  Also  included  on  each  figure  is  a  ’’thin  sheath  limit.” 
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Figure  2.  Separation  Curves  for  a  =  Figure  3.  Separation  Curves  for  a  = 

0.  01  and  Various  Values  of  N.  The  1.  0  and  Various  Values  of  N,  The 

dashed  curve  is  the  thin  sheath  limit  dashed  curve  is  the  thin  sheath  limit 


Figure  4.  Separation  Curves  for 
a  =  100  and  Various  Values  of  N. 
The  dashed  curve  is  the  thin  sheath 


This  is  obtained  by  plotting  0,  as  taken  from  the  a  =  0  separation  curve,  as  a 
function  of  0*  *  \P%-2 a.  The  reasoning  is  that  in  this  limit  the  particle  loses 

energy,  <*,  in  traversing  an  infinitesimally  thin  sheath  and  then  moves  in  a  field - 
free  region  starting  from  r  =  1,  with  an  initial  velocity  characteristic  of  the 
reduced  energy. 

One  might  expect  that:  (A)  separation  curves  for  successively  higher  powers 
of  N  would  form  a  monotonic  sequence,  and  (B)the  thin  sheath  limit  would  be 
approached  as  N  In  fact,  inspection  of  the  figures  shows  that  (A)  is  definitely 

not  true  while  (B)  may  or  may  not  be.  For  a  -  100,  the  curve  that  deviates  most 
from  the  thin  sheath  limit  (returns  most  particles  to  the  3phere)  is  that  for  N  =  2, 
while  for  a  =  1.0,  the  curve  for  N  =  3  plays  this  role.  Separation  curves  for 
ot  ~  0.  1  show  a  maximum  deviation  from  the  thin  sheath  limit  for  N  -  7.  For 
a  -  0.01,  the  separation  curves  for  values  of  N  up  to  12  form  a  monotonic  sequence, 
but  one  that  is  diverging  away  from  rather  than  converging  to  the  thin  sheath  limit. 

It  is  possible  that,  as  for  the  cases  with  larger  a,  there  will  be  a  reversal,  and 
that  for  some  higher  N  the  curves  will  start  to  approach  the  thin  sheath  limit;  but 
this  is  not  known. 

In  general,  changes  in  the  shape  of  the  separation  curves  corresponding  to 
changes  in  the  form  of  the  potential  seem  to  exhibit  some  complexity,  with  curves 
crossing  each  other  and  even,  on  occasion,  crossing  the  thin  sheath  limit.  Yet, 
in  spite  of  this,  a  general  statement  can  be  made:  unless  angles  of  emission  are 

high,  particles  having  energies  above  the  sphere  potential  will  have  a  good  chance 

2 

of  escaping.  For  example,  if  0  =4 a  (initial  kinetic  energy  =  twice  the  sphere 

potential)  and  N  =  12,  for  a  -  0.01  all  particles  having  9  <  1.05  radians  escape, 
for  a  -  1.0  all  particles  having  9  <:  1.40  radians  escape,  and  for  a  -  100,  par¬ 
ticles  having  9  <  1.  55  radians  escape. 


